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Abstract 

The experimental technique of fluorescence light observation is used in current 
and planned air shower experiments that aim at understanding the origin of ultra- 
high energy cosmic rays. In the fluorescence technique, the geometry of the shower 
is reconstructed from the correlation between arrival time and incident angle of 
the signals detected by the telescope. The calculation of the expected light arrival 
time used so far in shower reconstruction codes is based on several assumptions. 
Particularly, it is assumed that fluorescence photons are produced instantaneously 
during the passage of the shower front and that the fluorescence photons propagate 
on a straight line with vacuum speed of light towards the telescope. We investigate 
the validity of these assumptions, how to correct them, and the impact on recon- 
struction parameters when adopting realistic conditions. Depending on the relative 
orientation of the shower to the telescope, corrections can reach 100 ns in expected 
light arrival time, 0.1° in arrival direction and 5 g cm -2 in depth of shower maxi- 
mum. The findings are relevant also for the case of "hybrid" observations where the 
shower is registered simultaneously by fluorescence and surface detectors. 



1 Introduction 

Understanding the origin and nature of ultra-high energy (UHE) cosmic rays 
above 10 19 eV is a major challenge of astroparticle physics [1]. These cosmic 
rays are studied by detecting the atmospheric showers they initiate. Current 
and planned air shower experiments [2,3,4,5,6] use the technique of fluores- 
cence light observation: shower particles deposit energy in the atmosphere 
through ionisational energy loss. Part of this energy (of order 10~ 4 ) is emitted 
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isotropically at near-UV wavelengths in de-excitation processes. These fluo- 
rescence photons can be detected by appropriate telescope systems operating 
in clear nights. Typically pixel cameras with 25—100 ns timing resolution are 
used, where an individual pixel covers a field of view of about 1 — 1.5° in di- 
ameter (see e.g. Ref. [2]). The signal (light flux per time) is registered as a 
function of the viewing direction of the pixels. 

The first step to reconstruct the primary parameters of an observed air shower 
is given by the determination of the shower geometry. An accurate geometry 
reconstruction is, for instance, decisive for directional source searches; but it is 
also a prerequisite for reconstructing other important shower parameters such 
as the primary energy or the depth of shower maximum. We note that also 
the shower energies obtained from Auger ground array data are calibrated by 
the fluorescence telescopes [7]. 

The determination of the shower geometry is commonly performed in two steps 
in the fluorescence technique [9]. First, the "shower-detector-plane" (SDP) is 
determined as the plane spanned by the (signal- weighted) viewing directions of 
the triggered camera pixels (Fig. 1). Next, the geometry of the shower within 
this SDP is reconstructed based on the correlation between arrival time of the 
signals and viewing angle of the pixels projected into the SDP. The measured 
time-angle correlation is compared to the one expected for different shower 
geometries, and the best-fit geometry is determined. For the calculation of the 
expected time-angle correlation, the following function is in use (following e.g. 



where is the arrival time of the photons at camera pixel i (usually, a signal- 
weighted average arrival time is taken from the time sequence observed in a 
pixel), to is the time at which the shower axis vector passes the closest point 
to the telescope at a distance R p , c vac is the vacuum speed of light, Xo is the 
angle of incidence of the shower axis within the SDP, and Xi is the viewing 
angle of pixel % within the SDP (see also Fig. 1). Comparing the expected 
t~ Xi correlation to the observed one (i = l...n for n triggered pixels), the 
best-fit parameters R p , to and xo i n Eq. (1) are found by a ^-minimization. 
Together with the SDP derived previously, the shower geometry is then fully 
determined and can also be expressed in terms of shower impact point, arrival 
direction, and ground impact time. 

Eq. (1) is derived as follows. Assuming the fluorescence light to be emitted 
by a point-like object moving at c vac along the shower axis vector, the shower 
propagation time r s ho W er,i from point Si to the point at reference time t on 
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Fig. 1. Sketch of the shower geometry and quantities used in the derivations, 
the shower axis (Fig. 1) can be expressed as 

^shower, i ] 7 \~ • (^) 

c vac • tan(xo - Xi) 

Next, assuming the fluorescence photons to propagate on straight lines with 
c vac , the light propagation time r]i g ht,i from Si to the telescope is 

Tlight.i = r-f- r • (3) 

c vac • sm(xo - Xi) 

With Eqs. (2) and (3), and assuming an instantaneous emission of the fluo- 
rescence photons at Si, the expected arrival time U (relative to the time to of 
closest approach of the shower to the telescope) of fluorescence photons at a 
pixel viewing at an angle Xi becomes 



H ^0 ^shower, i flight, i 
Rn ( I 

— to + 



Cvac \ sin (Xo - Xi) tan(xc 
R p + fXo-Xi 



to + ^tan ^-^) (4) 

Cvac V * 



which equals Eq. (1). 

Thus, the derivation of Eq. (1) for calculating the expected time-angle corre- 
lation is based on the following assumptions: 
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• the spatial structure and the propagation of the shower disk can be approx- 
imated by a point-like object moving at c vac , 

• the fluorescence light is produced instantaneously, 

• the fluorescence light propagates with c vac , 

• the fluorescence light propagates on a straight line. 

In this article, we investigate the validity of these assumptions. The impact 
of the corrections on reconstruction parameters is studied. The results are 
relevant both for observations with fluorescence telescopes alone and for "hy- 
brid" observations where the shower is registered by fluorescence and surface 
detectors. 

The plan of the paper is as follows. In Section 2, we discuss the various assump- 
tions and their corrections individually. In Section 3, the impact on shower 
reconstruction is studied. Conclusions are given in Section 4. 



2 Analysis of individual effects 

We discuss step-by-step the individual effects given by 

• the spatial structure and speed of the shower disk (instead of a point-like 
object moving with c vac ), 

• the delayed (instead of instantaneous) fluorescence light emission, 

• the reduced propagation speed of light (instead of c vac ), 

• the bending of light (instead of straight-line propagation). 

2. 1 Spatial structure and speed of shower disk 

To check the assumption of the shower propagating as a point-like object 
with c vac on a straight line, one may first regard the fastest particles during 
the cascading process. Assuming, as a rough estimate, an elasticity of 50% 
per interaction, the energy of the leading particle in a hadronic air shower is 
E\ p ~ (Eq/Aq) -0.5 n after n interactions for a primary particle of energy Eq and 
mass Aq. For n = X max / A ~ 10 (the depth of shower maximum in units of the 
hadronic interaction length), the energy of the leading particle is ~10~ 3 i?o for 
primary protons and of order ~10~ 5 -E for primary iron. Hence, E\ p > 10 13 eV 
for primary energies E > 10 18 eV around shower maximum, which is the most 
relevant portion of the shower development for fluorescence light observations. 
In this case, the accumulated time delay of the leading particles with respect 
to an imaginary shower front moving with c vac from the first interaction to 
X miB is <^Cl ns. This is negligible compared to current timing resolutions of 
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giant shower detectors. Lateral deflections of these particles due to transverse 
momenta in interactions or deflection in the Earth's magnetic field are also 
sufficiently small (below ~1 m)Q For the case of UHE shower observations 
by fluorescence telescopes we conclude that the fastest shower particles can 
in reasonable approximation be assumed to move on a straight line along the 
shower axis with c vac - 

The main contribution to the fluorescence signal in the shower, however, is 
due to lower-energy secondaries, particularly electrons and positrons between 
0.1 MeV and several 100 MeV [11] These have larger lateral displacements 
from the shower axis and larger longitudinal time delays with respect to the 
shower front. 

Concerning the lateral width of the fluorescence shower beam, about 80% 
of the total fluorescence signal is produced within ~75 m around the shower 
axis [11]. The impact of the finite shower width on the fluorescence reconstruc- 
tion and how to correct it, was previously studied in detail [13]. It was shown in 
Ref. [13] that choosing too small a photon collection angle around the shower 
axis during reconstruction can lead to a signal loss and underestimation of the 
primary energy in nearby showers. 

Here we study the longitudinal time delay of secondaries using the CORSIKA 
code [14]. In Fig. 2 the time delay of secondaries, weighted according to their 
contribution to the energy release into air and thus to the fluorescence sig- 
nal, after the arrival time of the first particles is shown (10 19 eV shower at 
maximum, for particles closer than 75 m from the axis; results are practically 
identical for primary proton and iron showers). One can note a sharp ini- 
tial increase of the cumulative distribution (about 50% of energy is deposited 
within the first 3—4 ns after the fastest particle), with a long tail towards 
larger delays. The average time delay is ~8 ns, corresponding to a shower 
"thickness" of a few meters, which is in reasonable agreement with measure- 
ments of particle delays in air showers (see e.g. Ref. [15]). As can also be seen 
in Fig. 2, the distribution of time delays changes only marginally with the 
shower development stage. 

The delay of secondaries with respect to the fastest shower particles results in a 
small constant time offset of the observed shower compared to the assumption 



Time delay and lateral deflection of the leading particles may become non- 
negligible in case of considerably smaller Eq or larger n (the latter being rather 
relevant for ground array observations of near-horizontal showers). 
2 Note that for the energy transfer from > 0.1 MeV electrons to fluorescence pho- 
tons, the production of even lower-energy (e.g. < 1 keV) electrons is important (for 
instance, the cross-section for exciting the main molecular bands (cf. Section 2.2) 
has a sharp peak at about 20 eV electron energy). However, the additional delay 
from this intermediate step is ^ 1 ns and, thus, negligible for this analysis [12]. 



5 




time delay [ns] 

Fig. 2. Cumulative energy deposit (normalized to unity) as a function of time delay 
with respect to the fastest particle. The plot refers to a 10 19 eV proton at shower 
maximum (and at 100 g cm -2 smaller /larger depths) and includes particles within 
75 m of the shower axis. The average time delay is indicated be the dotted line. 
Simulations were performed with CORSIKA [14] / QGSJET 01 [16]. 

of the shower moving with c vac . This might be less relevant for observations 
with fluorescence telescopes alone, since in this case, only the relative timing 
between the pixels is used to determine the spatial shower geometry. For hy- 
brid observations, however, usually the arrival time of the first particle in the 
ground detector is taken, while in fluorescence telescopes, usually an average 
time from a fit to the signal viewed by a pixel is used. Then, comparing the tim- 
ing signals from ground and fluorescence detectors, the small shift due to the 
finite shower thickness should be taken into account The precise value of the 
delay will depend on the specific procedure of signal extraction applied during 
reconstruction. As a rough estimate, the delay is of order r t hi C k — 5—6 ns. 

To summarize, the leading particles in >10 18 eV showers can be considered to 
propagate along the shower axis with c vac , and one can set Ti ea di ng ,i — T s hower,i 
with r shower j given by Eq. (2). Compared to these particles, the secondaries 
relevant for the fluorescence light are slightly delayed due to the finite shower 
thickness by r thick , i.e. this term has to be added on the r.h.s. of Eq. (4). 

2.2 Fluorescence light production 

During propagation, the shower particles excite and ionize air molecules. Fluo- 
rescence light is then emitted by de-excitation and recombination. Most of the 

3 For ground detectors located at larger distances from the shower axis, the curva- 
ture of the shower front needs to be accounted for in addition. 
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fluorescence light originates from transitions from the second positive system 
(2P) of molecular nitrogen N 2 and the first negative system (IN) of ionized 
nitrogen molecules [8]. 

Typical excitation times are of the order ~10 -6 ns [17] and negligible for cur- 
rent fluorescence telescopes. De-excitation times, in turn, can exceed 30 ns. 
Depending on the local atmospheric conditions and on the specific transition 
system, quenching processes (radiationless transitions by collisions with other 
molecules) can substantially reduce the mean de-excitation time of the radia- 
tive processes. 

The total reciprocal lifetime 1 / r v > (p, T) of an electronic vibrational state v' 
can be expressed as a function of pressure p and temperature T as (see e.g. 
Ref. [18] and references therein) 

-J— = J-(i + _2_ 

v(p,T) t % , \ p' v ,(T) 

Here, 1/tq , is the reciprocal intrinsic lifetime defined as the sum of all con- 
stant transition probabilities and p' v , (T) is a reference pressure for a given gas 
mixture defined as the pressure where the collisional deactivation constant 
equals the reciprocal intrinsic lifetime [18]. 

Fig. 3 shows the calculated lifetimes as a function of height for the three main 
sets of bands 2P(0,z/'), 2P(l,z/') and lN(0,z/'), assuming dry air (78.1% N 2 , 
20.9% 2 and 1% Ar) and temperature profiles corresponding to the typical 
conditions at the Auger Observatory [19]. Also shown is the averaged lifetime, 
weighting the emission bands according to their relative (altitude dependent) 
intensities. The width of the weighted line indicates the effect of an arbitrary 
temperature variation of ±40 K to show the minor dependence of the averaged 
lifetime on reasonable variations of the actual atmospheric conditions. At very 
high altitudes of 30—40 km, the averaged lifetime is 15—25 ns. With decreasing 
altitude, the quenching effect reduces the lifetime; thus, in general, the delay 
of fluorescence light emission with respect to the passing shower front is a 
differential effect that changes during the shower development (smaller delay 
deeper in the atmosphere) o At heights below ~20 km where showers are 
typically observed by ground-based observatories, lifetimes of a few ns are 
reached. 

The average lifetime Td ee xc [i n ns] (weighted line in Fig. 3) can in good approx- 
imation be parameterized as a function of height h a.s.l. [in m] of the emission 



4 Anecdotally, this means the front of fluorescence light emission can move with an 
apparent velocity larger than c vac through the atmosphere. 
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Fig. 3. Lifetime of the three main sets of bands as a function of height a.s.l. for dry 
air. The thick line shows the averaged lifetime, weighted according to different inten- 
sity fractions. The width of the line indicates the effect of a change in temperature 
by ±40 K. 

point by 

TdeexcO) = . : , (6) 

a ■ e n i n + 1 



with r = 37.5 ns, H = 8005 m and a = 95. The term Td e0X c(^) has to be 
added to the r.h.s. of Eq. (4). 



2. 3 Reduced speed of light 



The propagation speed of light v = c va _ c /n is reduced compared to the vacuum 
case by the local index of refraction of air n. The change of n with wavelength 
is small (<3%) [20] within the fluorescence window of about 300—400 nm. 
Following Ref. [21], the index of refraction can be parametrized as a function 
of altitude h as 

n(h) = l + (n -l)^ (7) 
Po 



with the atmospheric density profile p(h); Uq and po are the reference values at 
sea level. The propagation time of refracted light over a small line element ds 
is then given by dr Te{r ~ n(h)ds/c vac and for propagation between two points 
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distance [km] 

Fig. 4. Arrival time difference (^reai — ^vacuum 

) due to the effect of reduced speed of 
light. The telescope is placed at 1.4 km a.s.l. corresponding to the altitude of the 
Auger telescope station "Los Leones." 



at altitudes h 2 and hi (h 2 > hi) by 



1 h2 

Trefr = q / 7l(h) dfl 

C vac COS V J 
hi 



with i? being the zenith angle of the propagation direction of light. 

In Fig. 4, the difference of light arrival times (between the cases of vacuum 
and reduced speed of light) is shown as a function of the location of emission 
point with respect to a telescope. The parametrization of p(h) is taken from 
Ref. [19] for the example of the southern Auger Observatory. As expected, 
for fixed distance between emission point and telescope, time differences grow 
for propagation closer to ground due to the larger value of n. Differences of 
20—25 ns or more can occur. For a single air shower, the effect changes along 
the longitudinal shower path, depending also on the relative orientation of 
shower axis and telescope. For instance, the time difference along the shower 
path typically varies less for showers pointing towards the telescope. 

In Eq. (4), Ti ight) j is replaced by r re f r)i . For convenience, one can express r rc f r j 
using Eq. (3) by replacing c vac with c re f r ,i = s/r re f r i , defined as the effective 
speed of refracted light along the path of length s between emission point and 
telescope. 
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2.4 Bending of light 



Due to refraction, the emitted light propagates on a bent trajectory. In turn, 
the direction of the incidence angle of the observed light does not point towards 
the real emission point, see Fig. 5 (right). More specifically, the zenith angle 
of down-going light is continuously reduced during propagation!!] The zenith 
angle difference A$ = $ rea i — ^a PP > between the observed light direction 
^ap P (towards the apparent emission point) and the straight-line direction i9 rea i 
(towards the real emission point) has been calculated from ray tracing with 
n(h) from Eq. (7); it is shown in Fig. 5 (left) as a function of the position of 
the emission point in the atmosphere relative to the telescope. As an example, 
an angular difference of 0.02° implies a ~12 m upward shift of the apparent 
emission point for a vertical shower at 30 km distance which corresponds 
to a ~40 ns shift in time. These shifts change over the longitudinal viewing 
direction towards an air shower. In case of hybrid observations where timing 
signals of fluorescence and ground detectors are combined, the impact time on 
ground estimated from the telescopes will be delayed compared to the actual 
one. 

For a vertical shower, or, more generally, for showers with ^sdp = 90° (cf. 
Fig. 1), Xi i n Eq. (4) is just reduced by Ai?, as the refracted light direction 
still points towards the actual shower axis. In general, however, this effect 
slightly shifts the refracted light signals out of the actual SDP, and this shift 
usually changes along the shower path. Thus, the apparent SDP (which, in fact 
need not be a "plane" anymore) may slightly be tilted compared to the real 
one. To still permit the practical approach of fitting the best shower geometry 
within a plane only (instead of testing the full phase space), the projected 
shift Ai9 • sin^sDP is taken as a correction. Thus, in Eq. (4), x% is replaced by 
X re fr,; — Xi — -sin t? SDP where Xrefr.i denotes the effective viewing direction of 
pixel i due to refraction. To account for the possible slight tilt of the apparent 
SDP, which is expected to be no larger than Ai9 max ~ (few times) 0.01°, the 
best-fit SDP might be found in an iterative procedure. 

Finally, we note that the additional time delay due to the increased, bent path 
length compared to the straight-line connection (see sketch in Fig. 5) is <^Cl ns 

5 We consider here only the case of a stable atmosphere with a standard decrease of 
p(h) and n{h) with height as given by Eq. (7) and Ref. [19]. We note, however, that 
the path of refracted light can become more complicated for specific atmospheric 
conditions such as atmospheric inversion, or in case of a strongly radiating ground 
leading to a local heating of air. The impact of the latter on the fluorescence tech- 
nique might be reduced due to the fact that observations are only performed well 
(~1— 2 h) after / before sunset; also, the shower path very close to ground usually 
is below the field of view of the telescope (~1° elevation of lower edge of field of 
view) . 
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distance [km] observer distance 



Fig. 5. Zenith angle difference A$ = # rea j — t9 app between direct and curved path 
due to light refraction as a function of the location of the emission point relative to 
the telescope. The telescope is placed at 1.4 km a.s.l. corresponding to the altitude 
of the Auger telescope station "Los Leones." 

and can thus be neglected. 



3 Impact on shower reconstruction 

Taking the discussed effects into account, Eq. (1) is finally replaced by 



U — to —~ / r H r + Tthick + Tdeoxc.j (9) 

Cvac tan(Xo - Xrefr.j) C reir ^ Sin(x - Xrefr,i) 

The index i indicates that these quantities, for a given shower geometry, de- 
pend on the viewing direction of pixel i. One caveat, as discussed in Sec. 2.4, is 
that the bending of light slightly changes the apparent SDP (within which the 
angles xo an d Xrefr,; are defined). It is worthwhile to note that all correction 
terms depend only on shower geometry but not on shower physics such as 
the primary particle type, which facilitates their application in shower recon- 
struction codes. r thick can, to a good degree, be treated as a constant; T deexCji 
depends on the altitude of the emission point; and c re f r) j and Xrefr.i depend on 
the locations of emission point and telescope. 

The time shifts introduced by the various effects along the viewing direction \i 
towards the shower are displayed in Fig. 6 for different shower geometries. The 
distance between impact point and telescope were fixed to 15 km (thin line) 
and 40 km (thick line), and for each distance three different shower inclinations 
of xo = 50°, 90°, 130° are considered. Here, for simplicity $sdp = 90° is taken 
such that 1 90° — \o\ is identical to the shower zenith angle. In this case, the 
effect from light bending is minimized concerning the change of the SDP and 
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maximized concerning Xi — Xrefr,i- 

In Fig. 6 (a), the overall shapes of t; L vs. x% are given, which differ for the 
different geometries. The shift of the arrival times, compared to the previous 
approach, is shown in Fig. 6 (b) when taking all effects into account. The 
contributions from the individual effects are provided in Figs. 6 (c)— (e). For 
the bending of light, in Fig. 6 (f), also the shift between apparent and effective 
viewing angle is given. One sees that the time delays are geometry dependent 
and can reach, and even exceed, 50—100 ns. 

One also sees in Fig. 6 that the time delays change along the shower track 
in an individual event. When reconstructing the shower as a whole, the fit- 
ting procedure then minimizes the overall x 2 by adjusting simultaneously R p , 
t and Xo- To investigate the effective impact of the corrections on the final 
reconstruction parameters, events were generated using CORSIKA [14] with 
the hadronic interaction model QGSJET 01 [16]. The shower sample consists 
of proton induced showers with energies of 10 18 , 10 19 and 10 20 eV and zenith 
angles of 0, 45 and 60 deg (100 events per combination with random azimuth 
angles). The detector simulation and the event reconstruction was performed 
using the Auger software package described in [22,23]. These data were re- 
constructed with and without accounting for the discussed effects; or, more 
specifically, using once Eq. (1) and once Eq. (9) in the reconstruction, and 
comparing the differences. 

As the time delays from individual parts of the shower were found to depend 
on distance and relative orientation of the shower to the telescope, we plot 
in Fig. 7 the change in the parameters R p , to and xo as a function of the 
minimum viewing angle (MVA), divided in different distance bins. The MVA 
is defined as the smallest angle under which the reconstructed air shower is 
seen by the telescope. Some dependence on MVA and distance can indeed be 
seen, as expected also from projection effects (for instance, a given angular 
offset in viewing direction leads to a larger shift along the shower axis for 
showers pointing towards the telescope than for vertical showers), or from an 
accumulation of certain effects with distance (such as the time delay due to the 
reduced speed of light). The actual impact of the corrections on an individual 
event is more complex, however, and has some dependence also on parameters 
other than geometry. For instance, the shower track of a higher-energy event 
can be observed out to larger altitudes due to the increased light output, such 
that these parts of the shower track can also contribute to the geometry fit. 
Thus, an a posteriori correction of the geometry parameters determined with 
Eq. (1) is not straightforward and would neglect individual event properties. 

For values of the MVA larger than 40—50°, deviations up to about 15 m (in 
Rp), 30 ns (in t ), and 0.04° (in xo) are observed. Differences in R p are typically 
larger in case of more distant showers. At smaller values of the MVA (showers 
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bending of light 
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Fig. 6. Upper left plot: light arrival time U vs. light arrival angle (or pixel viewing 
direction) Xi f° r different shower geometries (thick black (thin red) lines: shower 
impact point at 40 km (15 km) distance from the telescope; shower inclination xo 
= 130° (solid), 90° (dashed), 50° (dotted); in all cases i9 S dp = 90°; shower track 
shown up to 50 km distance from the telescope). Upper right to lower left plot: 
differences in expected light arrival time for the given shower geometries between 
old and new reconstruction including all effects (upper right) and for individual 
effects gned. Lower right plot: differences of viewing angles towards apparent 

and actual emission point due to refraction. 

pointing towards the telescope), also larger deviations are possible, e.g. shifts 
in Xo of 0.1° or more. In terms of differences in arrival directions (the relevant 
quantity for directional source searches), differences are typically around 0.05°, 
but can exceed 0.1°. A systematic shift can be noted to slightly overestimate 
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Fig. 7. Scatter plot of differences between old and new reconstruction in the geom- 
etry parameters R p , to and xo as a function of the minimum viewing angle. The 
events are divided according to their distance d of the impact point to the telescope 
(black dot: d < 20 km, green dot: d = [20, 33] km and red open circle: (i > 33 km). 
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Fig. 8. Differences between old and new reconstruction in shower zenith angle as a 
function of the minimum viewing angle (dots indicate individual events, red stars 
the average value). 

the shower zenith angles when neglecting the discussed effects, see Fig. 8. Shifts 
in energy are usually small (~ 0.5—1% on average). Reconstructed values for 
the depth of shower maximum are typically changed by 2—3 g cm -2 , with a 
trend of the corrected X max values being increased, and with larger corrections 
(5 g cm -2 and more) towards smaller values of MVA. 



4 Conclusion 

The assumptions used in the "classical" function of Eq. (1) for reconstructing 
the shower geometry from fluorescence light observations were investigated. 
The finite shower thickness leads to an energy deposition in air by secondaries 
which is delayed, compared to the shower front, by about 5—6 ns (with some 
dependence on the specific light collection algorithm employed). The emission 
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of fluorescence light is further delayed due to the finite lifetime of the transi- 
tions which, due to quenching, is altitude dependent. Typical values are a few 
nanoseconds up to 20 km height, and >15 ns for heights above 30 km. The 
propagation speed of light is reduced by the index of refraction of air. The 
delay, compared to a propagation with vacuum speed of light, depends on the 
locations of emission point and telescope, and can exceed 20—25 ns. Finally, 
another effect of refraction is the bending of light, which also depends on the 
locations of emission point and telescope. Angular differences between the ap- 
parent and actual emission point of 0.02° can occur, which may correspond 
to time shifts of several 10 ns. This effect can also lead to a slight tilt of the 
SDR 

All these corrections can be considered as geometrical ones, i.e. they are inde- 
pendent of specific properties of the individual showers other than their geom- 
etry. The corrected function for geometry reconstruction is given in Eq. (9). 
Compared to the previous approach, which assumed maximum propagation 
speed of both light and particles as well as no other delays, the effects of 
delayed timing (including the effect of bending of light) accumulate. In to- 
tal, differences of up to ~100 ns in estimated light arrival time are possible. 
Air shower experiments with comparable, or better, time resolution should 
take these effects into account. This refers both to data reconstruction and 
to implementing these effects in the shower— detector simulation. In terms of 
overall shower reconstruction parameters, corrections are typically 0.03—0.05° 
in arrival direction (with a systematic trend of overestimating the zenith angle 
when neglecting the effect), ~ 0.5—1% in energy and 2—3 g cm~ 2 in X max , 
but may in some cases exceed 0.1° and 5 g cm~ 2 . This is to be compared 
to typical reconstruction accuracies of ~0.6° (directional resolution) [24] and 
~11 g cm -2 (systematic X max uncertainty) [25] in case of Auger hybrid events. 

The increase in computing time for event reconstruction is modest, partic- 
ularly when applying the corresponding corrections only when approaching 
convergence in the minimization process (increase of ~20% or less, depend- 
ing on implementation). Some of the effects investigated in this work might 
be relevant also for shower detection techniques other than fluorescence tele- 
scope observations at ultra-high energy, e.g. Cherenkov light observations of 
air showers. 
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